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Abstract 

Independent of assumptions about the form of the quark-quark scattering 
kernel, K, we derive the explicit relation between the flavour-nonsinglet pseu- 
doscalar meson Bethe-Salpeter amplitude, Th, and the dressed-quark prop- 
agator in the chiral limit. In addition to a term proportional to 75, Th 
necessarily contains qualitatively and quantitatively important terms propor- 
tional to 75 7 • P and 75 7 • k k ■ P, where P is the total momentum of the 
bound state. The axial-vector vertex contains a bound state pole described 
by Tfj, whose residue is the leptonic decay constant for the bound state. The 
pseudoscalar vertex also contains such a bound state pole and, in the chi- 
ral limit, the residue of this pole is related to the vacuum quark condensate. 
The axial- vector Ward-Takahashi identity relates these pole residues; with the 
Gell-Mann-Oakes-Renner relation a corollary of this identity. The dominant 
ultraviolet asymptotic behaviour of the scalar functions in the meson Bethe- 
Salpeter amplitude is fully determined by the behaviour of the chiral limit 
quark mass function, and is characteristic of the QCD renormalisation group. 
The rainbow-ladder Ansatz for K, with a simple model for the dressed-quark- 
quark interaction, is used to illustrate and elucidate these general results. The 
model preserves the one-loop renormalisation group structure of QCD. The 
numerical studies also provide a means of exploring procedures for solving the 
Bethe-Salpeter equation without a three-dimensional reduction. 

Pacs Numbers: 14.40.Aq, 24.85.+p, ll.10.St, 12.38.Lg 
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I. INTRODUCTION 



7T- and 7^-mesons are the lightest hadrons and hence they play a significant role in 
the phenomenology of low-to-intermediate energy nuclear physics as mediators of the long- 
range part of the hadron-hadron interaction. They are easily produced in electron-nucleon 
and nucleon-nucleon collisions and therefore provide an ideal means of exploring models of 
hadronic structure and subnucleonic degrees of freedom in nuclei. As mesons, the simple 
quark- ant iquark valence-quark content of the 7r and K makes them the simplest light-quark 
systems one can study as strong interaction bound states; and this is a necessary step in 
developing a detailed understanding of their properties and interactions in terms of the 
elementary degrees of freedom in QCD. 

Mesonic bound states are described by the homogeneous Bethe-Salpeter equation [BSE], 
which is one of the Dyson-Schwinger equations [DSEs] characterising QCD. The homo- 
geneous BSE is an eigenvalue problem whose eigenvalue is P 2 , the square of the bound 
state mass, and whose eigenvector is the bound state amplitude [fully amputated, quark- 
ant iquark-meson vertex]. This bound state, or Bethe-Salpeter, amplitude is a crucial element 
in the calculation of production and scattering processes involving mesons, as illustrated in 
Refs. [0-§]. The BSE is familiar in the study of scattering and binding in two-nucleon sys- 
tems; and it is often illustrated, and its features explored, via the problem of two elementary 
scalars interacting via the exchange of a different elementary scalar ||. There have been 
many applications to the strong interaction meson spectrum, with recent, extensive studies 
in this general framework being those of Refs. which also cite related research. 

Bethe-Salpeter equation studies can be characterised by their treatment of the quark- 
antiquark scattering kernel, K(q, k; P); a concrete calculation being specified by a trunca- 
tion of the skeleton expansion for K. This kernel also appears implicitly in the DSE for the 
dressed-quark propagator [the QCD "gap equation"] via the dressed-quark-gluon vertex, 
(<?,£>)• In studies of the spectrum and interactions of bound states of light-quarks, where 
dynamical chiral symmetry breaking [DCSB] and Goldstone's theorem are particularly im- 
portant, it is crucial to ensure that K and are "mutually consistent", by which we mean 
that they must be such as to guarantee the preservation of the axial- vector Ward-Takahashi 
identity 0. Otherwise, as discussed and exemplified in Refs. [|f[^|,|10|, a qualitatively correct 
description of the light-quark meson spectrum is not possible; i.e., "fine-tuning" is necessary 
to properly describe the theoretical ideal of the chiral limit, and the observational fact that 
the pion is so- much-lighter than the characteristic hadronic scale: m p /2 ~ mAr/3 = M q , the 
constituent quark mass, but m n /2 ~ 0.2M g . 

The rainbow-ladder truncation of the quark-DSE and meson-BSE, without a three- 
dimensional reduction, is a specification of and K that ensures the preservation of the 
axial-vector Ward-Takahashi identity. It is fully specified by an Ansatz for the dressed- 
quark-quark interaction and allows a qualitatively and quantitatively good description of 
flavour-nonsinglet pseudoscalar, vector and axial vector mesons without fine-tuning, even in 
very simple models [ffTj . As such it is a phenomenologically efficacious tool in this sector. 

The fact that it describes the flavour- singlet pseudoscalar and scalar mesons poorly is 
not often mentioned. However, this defect is not crucial now that its source has been iden- 



tified and understood [12,13 1 . Employing a straightforward and systematic procedure for 



extending the rainbow-ladder truncation, a procedure that preserves the axial- vector Ward- 
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Takahashi identity at every order, allows one to analyse the attractive and repulsive terms 
order-by-order beyond ladder-truncation in the BSE. One finds, for example, that in the 
flavour-nonsinglet pseudoscalar channel the repulsive terms are approximately cancelled by 
attractive terms that arise at the same order, which explains why the ladder truncation 
provides a good approximation in this channel. This is not the case in the scalar channel 
where higher order terms do not cancel in a like manner but lead to a net repulsive effect, 
whose magnitude cannot be estimated a priori but which entails that the ladder truncation 
provides a poor approximation. In the flavour- singlet pseudoscalar channel, timelike gluon 
exchange diagrams arise when one improves upon the ladder-truncation and these provide a 
plausible mechanism for splitting the flavour- singlet and flavour-nonsinglet mesons; i.e., for 



generating a significant 97-77' mass splitting ||14|| . These observations illustrate and empha- 



sise that the rainbow-ladder truncation can lead to qualitatively and quantitatively reliable 
conclusions if used judiciously. 

Our goal herein is to provide a concrete illustration of the general results of Ref. ||; i.e., 
of the importance, feasibility and essential consequences of preserving the axial-vector Ward- 
Takahashi identity in BSE studies of quark- ant iquark bound states, and the extension of 
these results to SU (Nf > 3). We employ a renormalisable DSE-model of QCD that preserves 
the one-loop renormalisation characteristics of the dressed-quark and -gluon propagators and 
the quark-gluon vertex. This allows an explicit demonstration of the renormalisation group 
flow of the vacuum quark condensate, for example, and renormalisation point independence 
of physical observables in this framework. We concentrate on the ir- and i-T-mesons since 
this subsystem has all the complexity necessary for a complete discussion of the features we 
wish to elucidate. 

In Sec. II we describe the DSE for the renormalised dressed-quark propagator; this prop- 
agator being a critical element in the construction of the kernel in the BSE for meson bound 
states. We discuss this BSE in Sec. Ill, along with the constraints entailed by preserv- 
ing the axial-vector Ward- Takahashi identity. In Sec. IV we report a model study of the 
quark DSE and meson BSE, and a range of 7T- and iT-meson observables, illustrating the 
model-independent results derived in the preceding sections. We summarise and conclude 
in Sec. V. 



II. QUARK DYSON-SCHWINGER EQUATION 

In a Euclidean space formulation, with {7^,7^} = 25^, 7^ = 7 M and a-b = Y^i=i a ii>i, 
the DSE for the renormalised dressed-quark propagator is 

Sip)' 1 = Z 2 {i 1 -p + m bm ) + Z x f g 2 D, u (p - q )^- llM S(q)T a u (q,p) , (1) 

where D^(k) is the renormalised dressed-gluon propagator, T"(g;p) is the renormalised 
dressed- quark-gluon vertex, my, m is the A-dependent current-quark bare mass that appears in 
the Lagrangian and J* = J A d 4 g/(27r) 4 represents mnemonically a translationally-invariant 
regularisation of the integral, with A the regularisation mass-scale. The final stage of any cal- 
culation is to remove the regularisation by taking the limit A — ► 00. The quark-gluon- vertex 
and quark wave function renormalisation constants, Zi(/i 2 ,A 2 ) and Z 2 (/i 2 ,A 2 ) respectively, 
depend on the renormalisation point and the regularisation mass-scale, as does the mass 
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renormalisation constant Z m (/i 2 ,A 2 ) = Z 2 (fi 2 , A 2 ) _1 Z 4 (/i 2 , A 2 ). In Eq. ([[]), S, and mb m 
depend on the quark flavour, although we have not indicated this explicitly. However, in our 
analysis we assume, and employ, a flavour- independent renormalisation scheme and hence 
all the renormalisation constants are flavour- independent. 



A. General remarks about renormalisation 



1. Dressed- quark propagator 
The solution of Eq. ([!]) has the general form 



S{p)- 1 = i 1 -pA(p 2 ^ 2 ) + B(p 2 ^ 2 ) 



1 



Z{P 2 ^ 2 ) 

renormalised such that at some largef] spacelike-/^ 2 



«7 • p + M{p 2 , /i 2 ) 



Sip)' 1 



ry • p + m(/i) 



(2) 



(3) 



where m(/i) is the renormalised quark mass at the scale [i. In the presence of an explicit, chi- 
ral symmetry breaking, current-quark mass one has Z±m(\i) = ^mbm, neglecting 0(1/^ 2 ) 
corrections associated with dynamical chiral symmetry breaking that are intrinsically non- 
perturbative in origin. 

Multiplicative renormalisability in QCD entails that 



A{p 2 ^ 2 ) Z 2 (p 2 ,A 2 ^ 



A(/i 2 ,/i 2 ) 



1 



(4) 



Such relations can be used as constraints on model studies of Eq. ([[]). Explicitly, at one-loop 
order in perturbation theory, 



Z 2 (^ 2 ,A 2 " 



a{A 2 
a(fi 2 



"IF 



(5) 



where 7^ = |£ and /3 X = Nf/3 — 11/2, with £ the gauge parameter and Nf the number of 
active quark flavours. At this order, 

a{Q 2 ) 



7T 



QCD 



(6) 



1 Herein, by "large" we mean fj, 2 very-much-greater than the renormalisation-group-invariant 
current-quark mass for the s-quark so as to ensure that, in our model calculations, the renor- 
malisation constants are flavour- independent to better than 1%. It is possible to employ a mod- 
ified subtraction scheme in which the renormalisation constants are exactly flavour- independent, 
however, it does not quantitatively affect our results and hence is an unnecessary complication |l5| . 
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Clearly, at one-loop in Landau gauge [£ = 0], A(p 2 ,fi 2 ) = 1, and a deviation from this 
result in a solution of Eq. (|l|) is a higher-loop effect. Such effects are always present in the 
self-consistent solution of Eq. ([!]). 

The ratio M(p 2 , /i 2 ) = B(p 2 , fi 2 )/A(p 2 , /i 2 ) is independent of the renormalisation point in 
perturbation theory; i.e., with /i/p, 



At one-loop order: 



M(p 2 ,/i 2 ) = M(p 2 ,/i 2 ) = M{p 2 ) , Vp 2 



(7) 



m(fj,) = M(> 2 



m 



1 \ 7n 



A 5 

"qcd 



where m is a renormalisation-point-independent current-quark mass and 7 r , 
is the anomalous dimension at this order; and 



Z m {fi 2 ,A 2 



a (A 2 
a(/x 2 



(8) 



12/(33-2iV7) 



(9) 



In QCD, 7 m is independent of the gauge parameter to all orders in perturbation theory and 
the chiral limit is defined by rh = 0. Dynamical chiral symmetry breaking is manifest when, 
for rh = 0, one obtains m(/i) ~ Ofd/yU 2 ) 7^ in solving Eq. (0), which is impossible at any 
finite order in perturbation theory.^] This is discussed and illustrated in Sec. [IV A| . 



2. Dressed- gluon propagator 

In a general covariant gauge the renormalised dressed-gluon propagator in Eq. (|l|) has 
the general form 

^w = (v-^) f! ^+e^, do) 

where d(k 2 ,n 2 ) = 1/[1 + U(k 2 , fi 2 )], with U(k 2 ,fi 2 ) the renormalised gluon vacuum polar- 
isation. The fact that the longitudinal [^-dependent] part of D^l^k) is not modified by 
interactions is the result of a Slavnov- Taylor identity in QCD: k^D^k) = £k u /k 2 . We 
note that Landau gauge is a fixed point of the renormalisation group; i.e., in Landau gauge 
the renormalisation-group-invariant gauge parameter is zero to all orders in perturbation 
theory; hence we employ this gauge in all numerical studies herein. 
Multiplicative renormalisability entails that 

d{k 2 ^ 2 ) _ Z 3 (fi 2 ,A 2 ) _ , 2 _ 1 
d(k 2 ,fi 2 ) Z 3 (/i 2 ,A 2 ) a[fX ^ } d{^ 2 ^ 2 )- 



2 The arguments presented herein cannot be applied in a straightforward fashion to models whose 
ultraviolet behaviour is that of quenched QED4, such as Ref. jlfJ], where the chiral limit can't be 
defined in this way. The difficulties encountered in such cases are illustrated in Ref. [17]. 
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At one-loop order in perturbation theory 



^s(/i 2 ,A 2 



a{A 2 ) 
a(/z 2 ) 



2L 

ft 



(12) 



where 71 = |iV> - ±(13 - 30- 



3. Dressed- quark- gluon vertex 
The renormalised dressed-quark-gluon vertex in Eq. (P is of the form 

T a u (k,p) = ^T 1/ (k J p). 



(13) 



As a fully amputated vertex, it is free of kinematic singularities. The general Lorentz 
structure of T v (k,p) is straightforward but lengthy, involving 12 distinct scalar form factors, 
and here we do not reproduce it fully: 



T v (k,p) = y v Fi(k,p,fi) + 



(14) 



but remark that Ref. [[DJ, pp. 80-83, and Refs. |I1||T9"| provide an elucidation of its structure, 
evaluation and properties. 

Renormalisability entails that only the form factor F\, associated with the 7^ tensor, 
is ultraviolet- divergent. By convention, and defining: /i(fc 2 ,/i 2 ) = F\(k, —k,fi) , T u (k,p) is 
renormalised such that at some large spacelike-/^ 2 



1 



Since the renormalisation is multiplicative, one has 



h(k\p?) Zi{p?, A 2 ) "V"-' /i(/i 2 ,// 2 ) ' 
At one-loop in perturbation theory the vertex renormalisation constant is 

~a(A 2 ) 



01 



(15) 



(16) 



(17) 



where 7 r = |[|(3 + + k). 



B. Model for the quark DSE 



In order to exemplify the results of Ref. || , which we reiterate and generalise in Sec. [Ill A| , 
we must know the form of D^ik) and T v {k,p), not only in the ultraviolet where pertur- 
bation theory is applicable, but also in the infrared, where perturbation theory fails and 
lattice simulations are affected by finite- volume artefacts. D^ u (k) and T v (k,p) satisfy DSEs. 
However, studies of these equations in QCD are rudimentary and are presently best used 
only to suggest qualitatively reliable Ansatze for these Schwinger functions. That is why all 
quantitative studies of the quark DSE to date have employed model forms of D^ v {k) and 
T u (k,p). 
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1. Abelian approximation 



To introduce one commonly used pair of Ansatze, we use Eqs. (|5[), fll^p and flT7| ) and 
observe that 



2tf 7i _ ^7r 
Pi Pi Pi 



1 . 



(18) 



Hence, on the kinematic domain for which Q 2 = (p — q) 2 ~ p 2 ~ q 2 is large and spacelike, the 
renormalised dressed-ladder kernel in the Bethe-Salpeter equation for the (fully-amputated) 
Bethe-Salpeter amplitude behaves as follows: 



g 2 (fi 2 ) D^(p- q) \Tl(p + ,q + )S (q + 
= 4vr a(Q 2 ) D^(p - q) [£ 7 „ S frcc (g + ) 



S(q-)ri(q-,p- 



(19) 



x 



S tTee (q. 



where P is the total quark-antiquark momentum, p + = p + rjpP and p_ = p — (1 — f]p)P [see 



Eq. (25)]. This observation, and the intimate relation between the kernel of the pseudoscalar 
BSE and the integrand in Eq. ([!]) [|T^], provides a means of understanding the origin of an 
often used Ansatz for D^ v (k); i.e., in Landau gauge, making the replacement 



g 2 D lxv (k) ^Ana{k 2 



[IV 



(k) 



in Eq. (]l|), and using the "rainbow approximation": 



(20) 



(21) 



The Ansatz expressed in Eq. (Effi ) is often described as the 'Abelian approximation" 
because the left- and right-hand-sides are equal in QED. In QCD, equality between the 
two sides of Eq. fl2"Dp cannot be obtained easily by a selective resummation of diagrams. 
As reviewed in Ref. |l|], Eqs. (5.1) to (5.8), it can only be achieved by enforcing equality 
between the renormalisation constants for the ghost-gluon vertex and ghost wave function: 
Z\ = Z%. 

A mutually consistent constraint, which follows from Z\ = Z% at a formal level, is to 
enforce the Abelian Ward identity Z\ = Z2. At one-loop this corresponds to neglecting the 
contribution of the 3-gluon vertex to T u , in which case 7r — > |£ = Jf- This additional 
constraint provides the basis for extensions of Eq. (£HJ); i.e., using Ansatze for T u that are 
consistent with the vector Ward-Takahashi identity in QED [p0|1 , such as Refs. |2T|-[23|]. 

The combination of Abelian and rainbow approximations [with Z\ — 1 = Z2] yields a 
mass function, M(p 2 ), with the "correct" one-loop anomalous dimension; i.e., 7 m in Eq. ([Bp 
in the case of explicit chiral symmetry breaking or (1 — j m ) in its absence ||24}| . However, 
other often used Ansatze for T u p!8|j25|] yield different and incorrect anomalous dimensions for 
M(p 2 ) |p6[ . This illustrates and emphasises that the anomalous dimension of the solution 
of Eq. ([!]) is sensitive to the details of the asymptotic behaviour of the Ansatze for the 
elements in the integrand. One role of the multiplicative renormalisation constant Z\ is to 
compensate for this. 
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2. Model for g 2 D llu (p-q)T l/ (q,p) 
Herein we employ a model for the kernel of Eq. (|I|) based on the Abelian approximation: 

Zi 9 2 D, v (p - q)^S(q)K(q,p) - £ Q{(p - qf) D^(p - g)y 7 ,%)y7, , (22) 
with the specification of the model complete once a form is chosen for the "effective coupling" 

One consideration underlying this Ansatz is that we wish to study subtractive renormal- 
isation in a DSE-model of QCD and it is not possible to determine Z\ without analysing 
the DSE for the dressed-quark-gluon vertex; a problem we postpone. Instead we explored 
various Ansatze for Y u and found that, with Q{k 2 ) = 4ira(k 2 ) for large-fc 2 , there was always 
at least one Ansatz for Z\ that led to the correct anomalous dimension for M(p 2 ). This 
interplay between the the renormalisation constant and the integral is manifest in QCD and 
Eq. ( p2[ ) is a simple means of implementing it. 

In choosing a form for Q{k 2 ) we noted that the behaviour of a(k 2 ) in the ultraviolet; 
i.e., for k 2 > 1-2 GeV 2 , is well described by perturbation theory. Constraints on the form of 
Q{k 2 ) in the infrared come from the DSE satisfied by the dressed-gluon propagator, D^{k). 
As summarised succinctly in Refs. [0,E^], qualitatively reliable studies of this equation 
indicate that the dressed-quark-quark interaction is significantly enhanced in the infrared 



such that on this domain it is well represented by an integrable singularity Pq| . Combining 
these observations with Eqs. (|T8|)-(|2"0"|), which illustrate the necessary interplay between the 
anomalous dimensions of each term in the integrand of Eq. ([!]) , motivates the Ansatz 

G{k2) ^D5 4 (k) + ^Dk 2 e~ k2 /- 2 + 4vr f —^f(k 2 ) , (23) 



ur i 



2 ln 



r + (i + P /A 2 QcD y 



with F(k 2 ) = [1 - exp(-k 2 /[4:m 2 })}/k 2 and r = e 2 - 1. (We use N f = 4 and A* r 



N f 

L QCD 

0.234 GeV in our numerical studies.) This is a simple modification of the form used in 
Ref. [0; one which preserves the one-loop renormalisation group behaviour of QCD in the 
quark DSE. 

The qualitative features of Eq. (|23|) are clear. The first term is an integrable infrared 
singularity |2Pj and the second is a finite-width approximation to S 4 (k), normalised such 
that it has the same J d 4 k as the first term. In this way we split the infrared singularity into 
the sum of a zero-width and a finite-width piece. The last term in Eq. (^) is proportional 
to a(k 2 )/k 2 at large spacelike-/c 2 and has no singularity on the real-fc 2 axis. 

There are ostensibly three parameters in Eq. (|23|) : D, uo and m t (D = 2^ m m 2 in Ref. [0). 
However, in our numerical studies, using a renormalisation point fi = 19 GeV, which is 
large enough to be in the perturbative domain, we fixed to = 0.3 GeV(= l/[.66fm]) and 
m t = 0.5GeV(= l/[.39fm]), and only varied D and the renormalised u/d- and s-current- 
quark masses in order to obtain a good description of low-energy 7r- and i^-meson properties. 
As shown below, this is achieved with 

D = 0.781 GeV 2 , m u/d (p) = 3.74 MeV , m s {n) = 82.5 MeV . (24) 
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(We do not consider isospin breaking effects herein.) We chose the quoted values of uo and 
m t primarily so as to ensure that Q{k 2 ) ~ 4iTa(k 2 ) for k 2 > 2 GeV 2 , as illustrated in Fig. |l|. 



This is sufficient for our present illustrative study, just as the form in Ref. [Rj] was sufficient 
therein. However, increasing sophistication and/or an exploration of a broader range of 
observables is likely to require a more careful treatment of this or other parametric forms. 

Evolved according to Eq. (|8|), postponing until Sec. |IV A| the discussion of whether this 
formula is appropriate, the "best fit" mass values in Eq. ( [2"4| ) correspond to rn^jf^ = 6.4 MeV 
and m] GeV = 140 MeV. The "best fit" values are sensitive to the behaviour of Q(k 2 ) for 
k 2 ~ 1-2 GeV 2 , which can be illustrated by their dependence on uo: increasing uo — > 1.5 uo, 
while maintaining a good fit to it- and i^-meson observables, requires a ~10% reduction 
in the value of these masses. With minor modifications of our parametrisation we can 
satisfy our phenomenological constraints using current-quark masses that are a factor of 



~ 1.5-2 smaller, as canvassed in Ref. [30]. We would have to apply tighter constraints in 
our phenomenological application to make a statement about the current-quark masses of 
light-quarks that is more accurate than this. These considerations do not affect the ratio 
of our fitted current-quark mass values, m s (/i) /m u /d{^) = 22.0, which is consistent with 
Refs. |6|,|n|] and the discussion of Ref. |S0| . 



III. PION AND KAON BETHE-SALPETER EQUATIONS 

The renormalised, homogeneous, pseudoscalar Bethe-Salpeter Equation (BSE) is 

[r^;P)L = [ A [XH(q;P)}srK r t °(q,k;P), (25) 

Jq 

where: H = n or K specifies the flavour-matrix structure of the amplitude; Xh(q] P) — 
S(q+)T H (q; P)S(q-), with S(q) = diag(S' u (g), S d (q), S s (q)); q + = q + rj P P, g_ = q - (1 - 
rjp) P, with P the total momentum of the bound state; and r,. . . ,u represent colour-, Dirac- 
and flavour-matrix indices. 

In Eq. (P5D, K^(q, k; P) is the renormalised, fully-amputated quark-antiquark scattering 
kernel, which also appears implicitly in Eq. ([]]) because it is the kernel in the inhomogeneous 
integral equation [DSE] satisfied by T u (q; p). K^(q, k; P) is a 4-point Schwinger function ob- 
tained as the sum of a countable infinity of skeleton diagrams. It is two-particle-irreducible, 
with respect to the quark-antiquark pair of lines and does not contain quark-antiquark to 
single gauge-boson annihilation diagrams, such as would describe the leptonic decay of a 
pseudoscalar meson.0 The complexity of K"[£(q, k; P) is one reason why quantitative stud- 
ies of the quark DSE currently employ Ansatze for D^k) and T u (k,p). As illustrated 
by Ref. 0, however, the complexity of K^(q, k; P) does not prevent one from analysing 
aspects of QCD in a model independent manner and proving general results that provide 
useful constraints on model studies of QCD. 



3 A connection between the fully-amputated quark-antiquark scattering amplitude: M = K + 
K(SS)K + . . . , and the Wilson loop is discussed in Ref. 
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Equation ( p5|) is an eigenvalue problem. Solutions exist only for particular, separated val- 
ues of P 2 ; and the eigenvector associated with each eigenvalue, the Bethe-Salpeter amplitude 
[BSA]: T H {k; P), is the one-particle-irreducible, fully-amputated quark-meson vertex. In the 
flavour-octet channels the solutions with the lowest eigenvalues are the tt- and ii'-mesons.Q 
The solution of Eq. (E5|) has the general form |32| 



T H (k;P)=T H lb 



iE H (k;P)+ 7 -PF H {k;P) 



(26) 



+ 7 • k k ■ PG H {k; P) + ov k»P v H H {k; P) 



where for bound states of constituents with equal current-quark masses the scalar functions 
E, F, G and H are even under k ■ P — > —k ■ P and, for example, T K+ = ~ (A 4 + iX 5 ), 
with {A J , j = 1...8} the S77(3)-flavour Gell-Mann matrices. The requirement that the 
bound state contribution to the fully-amputated quark- ant iquark scattering amplitude: M = 
K + K(SS)K + . . . , have unit residue leads to the canonical normalisation condition for the 
BSA: 



2P„ 



tr 
tr 

k 



r H (q;-P)^^r H (q;P)S( q ^ 



+ 



(27) 



T H (q; -P)S(q + )T H (q;P) 



dS(q^ 



dP„ 



+ 



F)K rs (n k- P\ 

[x H (q; -P)] /y ' J [x H (k; P)\ ut , 



where T H (k, — P)* = C^Y H {—k, —P)C, with C = 7274, the charge conjugation matrix, and 
X t denoting the matrix transpose of X. 

In Eq. (p5l), En(k; P) 7^ acts as a "source" in the equations for Fnik] P), Gn(k; P) and 
Hn{k] P) so that, in general, these subleading Dirac components of IV(/c; P) are nonzero. 

A. Chiral symmetry 



In studies of flavour-octet pseudoscalar mesons a good understanding of chiral symmetry, 
and its explicit and dynamical breaking, is crucial. These features are expressed in the 
renormalised axial-vector Ward-Takahashi identity [AV-WTI] 



- iP»Tl(k; P) = S- 1 ^)^ + j^S- 1 ^) - M W If (*; P) - if (*; P) M w , (28) 



where: Mu\ = diag(m u (/x), iridic), m s (//)); the renormalised axial-vector vertex is given by 



rt(fc;P) 



tit 



+ / [xUK P )\sr K r t :(q,k;P) 



ta 



(29) 



*We do not consider the //-meson because of its mixing with the if, which cannot be described in 



ladder approximation [11,12|; the truncation of the BSE employed in Sec III B 
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with xfjg; P) = S(q, 



)T£(q;P)S(q- 



); and the renormalised pseudoscalar vertex by 



rf(*;P) 



tu 



T 



H 



75" 



+ 



tu 



xf(q;P)]„K™(q,k;P) 



(30) 



with x^ili P) = <^(?+)rf (<?; P)«S(g_). Multiplicative renormalisability ensures that no new 
renormalisation constants appear in Eqs. (p9|) and (|30| ) P,p3|. 

Any study whose goal is a unified understanding of the properties of flavour-octet pseu- 
doscalar mesons and other hadronic bound states must ensure the preservation of the AV- 
WTI, which correlates the axial-vector vertex, pseudoscalar vertex, bound state amplitudes 
and quark propagators. 



1. Chiral limit 



Equation ( |2"B| ) is valid for all values of the renormalisation-group-invariant current-quark 
masses, m particular for the chiral-limit when M {tl) Tf(k; P) = diag(0, 0, 0) = Tf(k; P) M (m) . 
In this case the AV-WTI is 



iP^k; P) = s-\k + ) l5 ^- + ^s-\k_) . 



(31) 



As a straightforward generalisation of Ref . || , it follows from Eqs. @ and fl3~Tl) that in the 
chiral limit the axial-vector vertex has the form 



T 



H 



^F R (k; P) + 7 • kk^G R {k- P) - <v k u H R (k; P) 



(32) 



P, 



+ fl(k;P)+ f HJ ^T H (k;P), 



where: F R , G R , H R and f are regular as P 2 -> 0; P^fg^fc; P) ~ 0(P 2 ); T H (k; P) is the 
pseudoscalar BSA in Eq. (0); and the residue of the pseudoscalar pole in the axial- vector 
vertex is the leptonic decay constant: 



/flPu 



2 , 2 



tr 



P H ) 757^(9 + )r ff (g;P)5(g_ 



(33) 



with the trace over colour, Dirac and flavour indices. In addition the chiral limit AV-WTI 
entails 



f H E H (k;0) = B(k 2 ) , 
F R (k;0) + 2f H F H (k;0)=A(k 2 ), 
G R (k;0) + 2f H G H (k;0)=2A'(k 2 ) 
H R (k; 0) + 2 f H H H (k; 0) = , 



(34) 
(35) 
(36) 
(37) 



where A(k 2 ) and B{k 2 ) are the solutions of Eq. ([[J) in the chiral limit. 

As remarked above, in perturbation theory, B(k 2 ) = in the chiral limit. The ap- 
pearance of a P(/c 2 )-nonzero solution of Eq. ([!]) in the chiral limit signals DCSB: one has 
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dynamically generated a momentum-dependent quark mass term in the absence of a seed- 
mass. Equations (|32]) and (|3"4T)-(pfl) show that when chiral symmetry is dynamically broken: 

1) the homogeneous, flavour- nonsinglet, pseudoscalar BSE has a massless, P 2 = 0, solution; 

2) the BSA for the massless bound state has a term proportional to 75 alone, with the 
momentum-dependence of En{k; 0) completely determined by that of the scalar part of the 
quark self energy, in addition to terms proportional to other pseudoscalar Dirac structures, 
Fh, Gh and Hh, that are nonzero in general; and 3) the axial-vector vertex, T,?(fc;P), is 
dominated by the pseudoscalar bound state pole for P 2 ~ 0. The converse is also true. 

The relationship, in the chiral limit, between the normalisation of the pseudoscalar BSA 
and fn has often been discussed; for example, Refs 



Consider that if one chooses to 
normalise Th such that Eh{0; 0) = -B(O), and defines the BSA so normalised as T H H (k; P), 
then the right-hand-side of fl2"T|), evaluated with Th — > T H H (k; P), is equal to 2P fl N H , where 
Nh is a dimensioned constant. Using Eqs. fl34])-(|37D it is clear that in the chiral limit 



N. 



H 



H 



(38) 



However, in model studies to date, this result is not obtained unless one assumes A{k 2 ) = 1. 
It follows that any kernel which leads, via ([I]), to A{k 2 ) = 1 must also yield, F H = = Gh = 
Hh, if it preserves the AV-WTI. In realistic model studies, where A(k 2 ) ^ 1, the difference 
between the values of Nh and fn is an artefact of neglecting Fh, Gh and Hh in (p6|) [l|]. 



2. Explicit chiral symmetry breaking 



Again as a straightforward generalisation of Ref. ||, in the presence of explicit chiral 
symmetry breaking the AV-WTI, Eq. (|2"%|), entails that both the axial- vector and the pseu- 
doscalar vertices have a pseudoscalar pole; i.e., 



T 



H 



5/A"'' * I 2 ^ 5 



I^Fr (k; P) + 7 • kk^Glik; P) - <v K Hi (k; P) 



(39) 



+ rt(k;P)+l 



H 



M 2 T H (k;P) 



H 



and 



r?(fc;P) = ^75 



iSg(k;P)+ 1 -PF R T (k;P)+ 1 -kk-Pgg(k;P) 

1 



(40) 



+ r HWT ^r H (k;P) 



with: £f , Fg, T*, G H R , Ql, Hi, Hi and % regular as P 2 - -m 2 H ; P^k; P) ~ 0(P 
and 



2V 



(41) 



where /# is given by Eq. (0), with massive quark propagators in this case, and the residue 
of the pole in the pseudoscalar vertex is 
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ir H = Z A £ |tr ^T H )\ 5 S(q + )T H (q;P)S(q^ . (42) 

The factor Z4 on the right-hand-side depends on the gauge parameter, the regularisation 
mass-scale and the renormalisation point. This dependence is exactly that required to ensure 
that: 1) r H is finite in the limit A — > 00; 2) r H is gauge-parameter independent; and 3) the 
renormalisation point dependence of th is just such as to ensure that the right-hand-side of 
Eq. d41|) is renormalisation point independent. This is obvious at one-loop order, especially 
in Landau- gauge where Z 2 = 1 and hence Z 4 = Z m . 

In the chiral limit, using Eqs. (H) and (ll)-©, Eq. §|) yields 



r° H = - W (qq)l, ~ (qq)l = Z 4 {fi 2 ,A 2 )N c f tr Dirac [S^{q)} , (43) 



where the superscript "0" denotes that the quantity is evaluated in the chiral limit and 
(qq) ^, as defined here, is the chiral limit vacuum quark condensate, which is renormalisation- 
point dependent but independent of the gauge parameter and the regularisation mass-scale. 
Equation ( flOf ) is the statement that the chiral limit residue of the bound state pole in the 
flavour-nonsinglet pseudoscalar vertex is (— (qq) ^) / f% ■ 
From Eqs. (f|l|) and fl43|) one obtains immediately 

flml = - [mM + m d (fx)} + O (mfj (44) 
f 2 K+ m 2 K+ = - [mM + m s (/i)] + O (m 2 q ) , (45) 

which exemplify what is commonly known as the Gell-Mann-Oakes-Renner relation. 



We emphasise that the primary result, Eq. (|4T|), of which Eqs. (|44|) and fl45|) are corol- 
laries, is valid independent of the magnitude of rh q . We can rewrite it in the form 

f 2 H m 2 H = -(qq)^M H (46) 

where we have introduced the notation 

-(qq)" = fHr H , (47) 



in order to highlight the fact that, for nonzero current-quark masses, Eq. (|41|) does not 
involve a difference of vacuum massive-quark condensates; a phenomenological assumption 
often employed.0 



5 In QCD, the integral defining th diverges logarithmically, like the trace of the chiral limit quark 
propagator (vacuum quark condensate) , which is the reason why the right-hand-side of Eq. (|4l| ) is 
independent of the renormalisation point. This is unlike the trace of the m^O quark propagator, 
which diverges quadratically. 
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B. Model BSE 



In order to exemplify the results of Ref. M, which we have reiterated and generalised 



in Sec. [Ill A| , we must have an explicit form for the kernel K^(q,k;P) in Eq. (|25|). The 



form must be such as to preserve the AV-WTI, Eq. (P8|) , which requires a truncation of 
the skeleton expansion for K^(q, k; P) that is consistent with Eq. ([22]); our Ansatz for the 
kernel of Eq. ([!]) . The "ladder truncation" fulfills this requirement |l2|Jll 



KIM k-P) = - Q{{k - qf) D*r(k -q) I 7m 



(48) 



tr 



in which case Eq. (^Sj) becomes 

,A 
"1 



T H (k; P) + f G((k - qf) D^ c (k - q)^S(q + )T H (q; P)<%_)^ 7l/ = 



(49) 



and the normalisation condition, Eq. ([27]), simplifies because the last term vanishes when 
Kl£(q, k; P) is independent of P M . 



IV. NUMERICAL RESULTS 
A. Solution of the quark DSE 
Using Eqs. ([I]) and (p2]) our model quark DSE is 



S(p, fi)- 1 = Z 2 t^ -p + Z 4 mQi) + £'(p, A) , (50) 
with the regularised quark self energy 

S'(P,A) = j q A mP-q) 2 )D^p-q)^S(q)^ } (51) 

where Q{k 2 ) is given in Eq. (|23|) . Equation fl50D is a pair of coupled integral equations for 
the functions A(p 2 ,fi 2 ) and B(p 2 ,fi 2 ) defined in Eq. (|]). 

In the case of explicit chiral symmetry breaking, rh ^ 0, the renormalisation boundary 
condition of Eq. (^) is straightforward to implement. Writing 

E'(p, A) = %1 ■ p (A'(p 2 , A 2 ) - l) + B\p 2 , A 2 ) , (52) 

Eq. @ entails 

Z 2 (fi 2 , A 2 ) = 2 - A\fi 2 , A 2 ) and m(/i) = Z 2 (fj 2 , A 2 ) m bm (A 2 ) + B\fi 2 , A 2 ) (53) 



and hence 



A(p 2 , /i 2 ) = 1 + A'(p 2 , A 2 ) - A V, A 2 ) , (54) 
B(p 2 , v 2 ) = m{n) + B'(p 2 , A 2 ) - 5 V, A 2 ) . (55) 
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From Sec. |II A 1|, having fixed the solutions at a single renormalisation point, fx, their 



form at another point, fx, is given by 

S-\p,n) = h- P A(p\fi>) + Btf,j?) = P£^ S-\p,») . (56) 

Z 2 {fx 2 ,A 2 ) 

[Recall that M(p 2 ) is independent of the renormalisation point.] This feature is manifest in 
our solutions. It means that, in evolving the renormalisation point to p, the "1" in Eq. (|54|) 
is replaced by Z 2 (fx 2 , A 2 )/Z 2 (fx 2 , A 2 ), and the "m(/i)" in Eq. Q5^) by m(p)\ i.e., the "seeds" 
in the integral equation evolve according to the QCD renormalisation group. 

As also remarked in Sec. [II A 1| , the chiral limit in QCD is unambiguously defined by 
rh = 0. In this case there is no perturbative contribution to the scalar piece of the quark self 
energy, B(p 2 ,fx 2 ), and, in fact, there is no scalar, mass-like divergence in the perturbative 
calculation of the self energy. It follows that Z 2 (fi 2 , A 2 )mt, m (A 2 ) = , VA and, from Eqs. ( |53|) 
and (|55D , that there is no subtraction in the equation for B(p 2 ,fi 2 ); i.e., Eq. (|55|) becomes 

B(p 2 ,fx 2 ) = B'(p 2 ,A 2 ) } (57) 

with liniA->oo B'(p 2 , A 2 ) < oo.[] In terms of the renormalised current-quark mass the existence 
of DCSB means that, in the chiral limit, M(fx 2 ) ~ 0(l//i 2 ), up to In /i 2 -corrections. 

In Fig. |2| we present the renormalised dressed-quark mass function, M(p 2 ), obtained by 
solving Eq. (|50|) using the parameters in Eq. (|24D , and in the chiral limit. [Recall that 
fx = 19 GeV, which is large enough to be in the perturbative domain.] 

It is clear from this figure that the light-quark mass function is characterised by a sig- 
nificant infrared enhancement, a direct result of that in the effective coupling, Q(k 2 ). In- 
troducing the Euclidean constituent-quark mass, M E , as the solution of p 2 = M 2 (p 2 ), the 
ratio: Mf /m/(/x), where / labels the quark flavour and mf(fx) is given in Eq. (p4|) , is a 
single, indicative and quantitative measure of the nonperturbative effects of gluon-dressing 
on the quark propagator. We find 

M E 

m f {fi) 

chiral limit 0.55 GeV oo (58) 
u/d 0.56 150 

s 0.70 8.5 

which clearly indicates the magnitude of this effect for light-quarks. The ratio M E /mf(fi) 
takes a value of 0(1) for heavy-quarks |37[] because the current-quark mass is much larger 



than the mass-scale characterising the infrared enhancement in the effective coupling, Aqcd- 
This means that in the spacelike region the momentum-dependence of the heavy-quark mass 
function is dominated by perturbative effects.^ 



This is a model-independent statement; i.e., it is true in any study that preserves at least the 
one-loop renormalisation group behaviour of QCD. 

7 Quark confinement entails that there is no "pole-mass" |35|, which would be the solution of 
p 2 + M 2 (p 2 ) = 0. Hence, this definition of M E is arbitrary; a factor of 2 is certainly unimportant 
with respect to the qualitative features that this quantity characterises. 
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In typical quark model calculations |36j the "constituent-quark" masses are M u / d = 
0.33 GeV and M s = 0.55 GeV. These are within a factor of 2 of the values in Eq. (p8|) ; 
and are in the ratio M u/d /M s = 0.60. From Eq.(§8[) we find M® /d /Mf = 0.80. The 
comparison of numerous DSE studies makes it clear that this correspondence between Mf 
and Mj? is robust. It provides a qualitative understanding of the nature of the "constituent- 
quark" mass; i.e., it is a quantitative measure of the nonperturbative modification of quark 
propagation characteristics by gluon dressing. Its magnitude is a signal of the enhancement 
of the quark-quark interaction in the infrared. 

The qualitative difference between the behaviour of M(p 2 ) in the chiral limit and in the 
presence of explicit chiral symmetry breaking is manifest in Fig. 0. In the presence of explicit 
chiral symmetry breaking Eq. (||) describes the form of M(p 2 ) for p 2 > 0(1 GeV 2 ). In the 
chiral limit, however, the ultraviolet behaviour is given by 



M(p 



2 x large_-p 2 



27T 2 7r , 



(qq)°) 



p z 



| In 



A 2 



1— 7n 



(59) 



where (qq)° is the renormalisation-point-independent vacuum quark condensate. f\ Analysing 
our chiral limit solution we find 



- (qq)° = (0.227 GeV) ; 



(60) 



This is a reliable means of determining (qq)° because corrections to Eq. (|59|) are suppressed 
by powers of A 2 QCB /fi 2 . 

Equation (E3) defines the renormalisation-point-dependent vacuum quark condensate 



(qq)l 



At=19 GeV 



= kl™o ^ 4 ( /i > A ) Nc J trDirac [ s rh=o(q)\ 



(0.275 GeV) 3 . (61) 



fi=19 GeV 



We have established explicitly that m(/i) (qq) ^ = constant, independent of /i with the value 
depending on the quark flavour, and hence 



(qq)°a - m (qqf 



(62) 



which unambiguously defines the renormalisation-point-independent current-quark masses. 
From this and Eqs. (p4|), (|60|) and (^1|) we extract the values of these masses appropriate to 
our model 



m u /d 



6.60 MeV, m, = 147 MeV. 



(63) 



Using Eq. (|8|) these values yield m u / d {jj) = 3.2 MeV and m s (fi) = 72 MeV, which are within 
~ 10% of our actual values in Eq. fl2H|). This indicates that higher-loop corrections to the 
one-loop formulae, which are present in the solution of the integral equation as made evident 



8 The momentum dependence of this result is characteristic of the QCD renormalisation group at 



one- loop [38 1 and demonstrates that the truncation we employ preserves this feature. 
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by A(p 2 ,fi 2 ) ^ 1, provide contributions of < 10% at p 2 = /i 2 . These contributions decrease 
with increasing p 2 ^\ 

From the renormalisation-point-invariant product in Eq. (|62"D we obtain 



- W ,=igcV = (ln [VAqcd]) 7 ™ {qq? = (0.241 GeV) 3 . (64) 

This result can be compared directly with the value of the quark condensate employed in 
contemporary phenomenological studies |39| : (0.236 ± 0.008 GeV) 3 . We note that increasing 
uj — > 1.5 u> in G(k 2 ) increases the calculated value in Eq. (|64T) by ~ 10%. Obtaining broad 
agreement with the contemporary phenomenological value of (9Q')° = iGeV was a means we 
employed to constrain the value of this parameter. However, we made no attempt to fine- 
tune uj nor thereby our calculated value of (<?g)°=iG e v- 

In conjunction with Eq. ( |64|) we define mj- GeV via Eq. (|) using Eq. (|63| ) 

miff = 5.5 MeV , m\ GeV = 130 MeV . (65) 



These values differ slightly from those discussed in Sec. |II B 2| for the reasons described above. 



It is now clear, from Eq. (|62|) , that lower values of the current-quark masses, as canvassed in 
Ref. |30|, are admissible in our phenomenological study only via an increase in (<??)°=iGeV- 



After this discussion of the vacuum quark condensate it is now straightforward to deter- 
mine the accuracy of Eqs. ( [44]) and ([EJ). Using experimental values on the left-hand-side, 
we find: 

(0.0924 x 0.1385) 2 = (0.113 GeV) 4 cf. (0.111 GeV) 4 = 2 x 0.0055 x 0.24 3 (66) 
(0.113 x 0.495) 2 = (0.237 GeV) 4 cf. (0.206 GeV) 4 = (0.0055 + 0.13) x 0.24 3 , (67) 

which indicates that 0(m 2 )-corrections begin to become important at current-quark masses 
near that of the s-quark. We emphasise that we did not use these equations in fitting the 
current-quark masses but, as described in Sec. |IV B| , solved the model BSE, Eq. ([49]), using 



the dressed-quark propagators obtained as a solution of Eq. (p0[) . In this procedure, changes 
in {qq)° effected by modifying the model parameters are compensated for by changes in m„/d 
and m s . Therefore the comparisons in Eqs. (^) and flS"?]), which involve a product of these 
quantities, remain meaningful and the conclusion remains valid.0 

B. Solution of the BSE 

To solve Eq. fl49|) we modify it by introducing an eigenvalue, A(P 2 ), 

T H (k; P) + A(P 2 ) J^g{{k - q) 2 ) Dfc{h - q)^S(q + )T H (q; P)S(g_)y7„ = , (68) 



9 Our model for the kernel of the quark-DSE is not constructed to preserve the two-loop, pertur- 
bative behaviour of QCD. Hence a direct comparison at this level is not meaningful. 

10 Using our calculated values of f n , m^, fx and m^, Tables | and |n[ the only change is in 
Eq. ©, where 0.237 -» 0.233. 
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which yields an equation that has a solution VP 2 , characterised by the value of A(P 2 ). The 
original problem is solved when that P 2 is found for which A(P 2 ) = 1, which will occur at 
P 2 < in our metric. 

For P 2 = — V 2 < 0, solving Eq. (|68|) requires the dressed-quark propagator <S(<?+) on 
the parabola: 4£ 2 P 2 3?(g 2 ) = S(g 2 ) 2 - 4(£ 2 P 2 ) 2 , in the complex-g 2 plane, and on the 

parabola: 4(1 — ^) 2 P 2 3?(g 2 ) = 3(g 2 ) 2 — 4((1 — £) 2 P 2 ) 2 , in the complex-g 2 plane. For complex 
arguments in the dressed-quark propagator the quark-DSE requires the effective coupling at 
complex values of its argument. Equation (^), however, still only requires Q(k 2 ) on the real- 
k 2 > axis. The specification of Q(k 2 ) is the primary element in the definition of the model 
and our Ansatz is motivated by studies that are restricted to real-fc 2 > 0. In employing 
this Ansatz at complex values of its arguments we are exploring an unconstrained domain. 
Solving numerically for S(p) in the complex-p 2 plane is straightforward. However, complex- 
conjugate branch points in the [confining] solution introduce numerical complications in 
solving the BSE for bound states containing a single heavy constituent 

The general form of the solution of the BSE, Eq. (f49"|), is given in Eq. (|26|) where the 
scalar functions depend on the variables k 2 and k ■ P, and are labelled by the eigenvalue 
P 2 . From this it is clear that the integrand in Eq. (^) depends on the scalars: k 2 , k ■ q, 
q 2 , q ■ P and P 2 , which takes a fixed-value at the solution; i.e., at each value of P 2 the 
kernel is a function of four, independent variables. Solving Eq. fl68|) can therefore require 
large-scale computing resources, especially since there are four, independent scalar functions 
in the general form of the solution. 

We employed two different techniques in solving Eq. fl68l) . In our primary procedure, 
(A), we treated the scalar functions directly as dependent on two, independent variables: 
E(k 2 , k ■ P; P 2 ), etc., which requires straightforward, multidimensional integration at every 
iteration. Storing the multidimensional kernel requires a large amount of computer memory 
but the iteration proceeds quickly. 

As an adjunct, and because we wish to elucidate qualitative and quantitative effects 
related to the k ■ P-dependence of the scalar functions in the BSA, we employed a Chebyshev 
decomposition procedure, (B). To implement this we write 

iVm ax 

E(k 2 ,k- P-P 2 ) « J2 'P(A; 2 ;P 2 )f/i(cos/3), (69) 

i=0 

with similar expansions for F, G = k ■ P G and H, where k-P = cos j3Vk 2 P 2 and {Ui(x); i = 
0, . . . , oo} are Chebyshev polynomials of the second kind orthonormalised according to: 

- f 1 dx Vl-x 2 Ui{x)UAx) = Su . (70) 

7T 7-1 
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This branch-point-pair is present because of the infrared enhancement in our Ansatz for the 



effective coupling but may be an artefact of the rainbow approximation, Eq. (|2l|). Ref. [22| demon- 
strates that dressing the vertex, consistent with the constraints of Refs. pTJ, can significantly affect 
the analytic properties of S(p) while maintaining the essence of quark confinement (35|; i.e., that 
S(p) not have a Lehmann representation. 
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Substituting this expansion into Eq. (|6q ) allows all-but-one of the integrals to be evaluated 
before beginning the iteration. One then solves for the Chebyshev moments, E(k 2 ; P 2 ). 
This procedure requires a large amount of time to set up the kernel but does not require 
large amounts of computer memory.^ 

In order to fit the parameter D in Q{k 2 ) and the current-quark mass m u /d{l^) we: 1) 
chose values for these parameters; 2) solved the quark DSE for S u /d(p)', 3) used Q{k 2 ) and 
the calculated form of S u m(p) to solve the pion BSE for the mass and BSA, T n (k; P), using 
procedure (A); and 4) used the calculated forms of S u /d(p) and r 7r (/c; P) to calculate f n from 
Eq. (P^D- We repeated this procedure until satisfactory values of m n and f n were obtained. 
Having thus fixed D we repeated the steps for the .ff -meson, varying m s {p) only, in order 
to obtain the best possible values of uik and fx- This led to the parameter values quoted 
in Eq. (B3) and the results listed in row one of Tables | and [[]J. This does not represent an 
exhaustive search of the available parameter space but is sufficient for our purposes. 

1. Discussion of the BSE solution 
From Tables |and||, and Eqs. (|24]), Q46|), (|47|) and (|64[), it is straightforward to calculate 



-<99>£=lGeV -(??)J=lGeV ~ (<?<?> ^=1 GeV 

(71) 

(0.245 GeV) 3 (0.284 GeV) 3 (0.317 GeV) 3 

showing that, for light pseudoscalars, the "in-meson condensate" that we have defined in- 
creases with increasing bound state mass; as does the leptonic decay constant, 
In Tables | and [TT] we list values of the dimensionless ratio 

n H = - w y/ 2 H . 72 

Jh™h 

A value of TZh = 1 means that Eq. (|4"T| ) is satisfied and hence so is the AV-WTI.Q Looking 
at the tabulated values of TZh it is clear that the scalar function H is not quantitatively 
important, with the AV-WTI being satisfied numerically with the retention of E, F and 
G in the pseudoscalar meson BSA. The values of 1Zh, and the other tabulated quantities, 



12 Equation ( |69|) is only an identity in the limit iV max — ► oo but, in the present example, an 
accurate representation of the solution is obtained with A^ max = 1 or 2, which is consistent with 
the observations of Ref . || . 

13 (—(qq)^)/fH is the residue of the bound state pole in the pseudoscalar vertex, just as fn is the 



residue of the bound state pole in the axial- vector vertex. As expected, {qq)Z = iQ e y ~ 



M =l GeV 



14 It illustrates that the pseudoscalar-meson pole in the axial-vector vertex is related to the 
pseudoscalar-meson pole in the pseudoscalar vertex in the manner we have elucidated. A finite 
value in the chiral limit emphasises that m 2 H oc A4h as Mr — > 0. 
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highlight the importance of F and G: F is the most important of these functions but G 
nevertheless provides a significant contribution, particularly for bound states of unequal- 
mass constituents. We note that a poor value of fn is tied to a poor value of IZh, which 
emphasises the importance of preserving the AV-WTI and hence Eq. We have checked 
explicitly that our complete solutions satisfy Eq. (j38|). 

The tables illustrate the rapid convergence of the Chebyshev decomposition procedure, 
(B), with accurate solutions being obtained with the zeroth and first moments for bound 
states of equal mass constituents.^ The same is true for bound states of unequal mass 
constituents provided all Dirac amplitudes are retained in the solution. 

One observes that when the k ■ P-dependence of the scalar functions in the meson Bethe- 
Salpeter amplitude is included, physical observables are independent of the momentum par- 
titioning parameter, rjp. These calculations elucidate the manner in which this necessary 
requirement in covariant bound state studies is realised: the bound state amplitude depends 
on rjp in just that fashion which ensures physical quantities do not. All the scalar functions 
in the BSA must be included to ensure this. 

The tables also illustrate clearly the effect of the current-quark mass via a comparison of 
the chiral limit pion with the physical pion, the kaon and a fictitious, pseudoscalar, ss bound 
state. One observes that fn is weakly sensitive to increasing the current-quark mass; for 
example, we obtain f° / ' f n = 0.97, which is consistent with expectations based on effective 
chiral Lagrangians. However, m 2 H rises quickly, before becoming sensitive to effects nonlinear 
in the current- quark mass.[^] 

We present the scalar functions in the BSA obtained as solutions of Eq. ( pE9"D in Figs. ||-|7|, 
focusing on the zeroth Chebyshev moment of each function, which is obtained via 

°E H (k 2 ) = - r dp sin 2 (3 U (cos (3) E H (k 2 , k ■ P; P 2 ) , (73) 

7T JO 

and similarly for F, G [G for the iC-meson] and H. Figure ^| illustrates that the momentum- 
space width of ^Euik 2 ) increases as the current-quark mass of the bound state constituents 
increases; Fig. f|, that °Fu(k 2 = 0) decreases with increasing current-quark mass but that 
°F H (k 2 ) is still larger at k 2 > 0.5 GeV 2 for bound states of higher mass; Fig. |], that °G H (k 2 ) 
[°GK{k 2 )} behaves similarly; and Fig. || that the same is true for .H#(fc;P) and that it 
is uniformly small in magnitude thereby explaining its quantitative insignificance. Where 
comparison is possible, these observations agree qualitatively with Refs. |||llf- 



15 We note that for equal-mass constituents G = k ■ P G is an odd function of k • P and hence 
°G = 0. Therefore it first contributes at 0(U\). 

16 We note that our calculated value of m S s is within 3% of that obtained for this fictitious bound 
state in Ref. Q. However, the results in Ref. Q depend on np, which is an artefact of the solution 
procedure adopted therein. In that study a derivative-expansion and extrapolation procedure was 
employed in order to avoid a direct solution for the quark propagator functions A(p 2 ) and B{p 2 ) 
at complex values of their arguments. In addition, in most instances only the zeroth Chebyshev 
moment was retained in the Chebyshev expansion of the scalar functions in the Bethe-Salpeter 
wave function: XH(q;P) = S(q + )T H (q; P)S(q~). 
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In Fig. |7| we illustrate the large-fc 2 behaviour of the scalar functions in the pseudoscalar 
BSA. The momentum dependence of Efj{k 2 ) at large- A; 2 is identical to that of the chiral- 
limit quark mass function, M(p 2 ), in Eq. ( |59| ) [ [40|j , and characterises the form of the quark- 
quark interaction in the ultraviolet. Figure elucidates that this is also true of °Fff(fc 2 ), 
k 2 °GH{k 2 ) [k 2 °GK(k 2 ) for the i^-meson] and k 2 if#(fc 2 ). Each of these functions reaches its 
ultraviolet limit by k 2 ~ 10 GeV 2 , which is very-much-less-than the renormalisation point, 
/x 2 = 361 GeV 2 . 



In order to verify Eqs. (|34])-(|3~7D it is necessary to consider the inhomogeneous axial- vector 



vertex equation, Eq. in our truncation of the DSEs, which is 

rf#; p) = z 2l5l ^ - £ Q((k - q f) D^(k - g)y7M«%+)rt(g; P)S(?-)y7* . (74) 

From the homogeneous BSE one already has the equations satisfied by E H (k]0), F H (k;0), 
G H (k;0) and H H (k;0). To proceed, one substitutes Eq. (0) for V^(k;P) in Eq. ([T§. 
Using the coupled equations for EH{k;0), etc., one can identify and eliminate each of the 
pole terms associated with pseudoscalar bound state. [It is the fact that the homogeneous 
BSE is linear in the BSA that allows this.] This yields a system of coupled equations for 
FR.(k;0), Gft(k;0) and H^(k;0), which can be solved without complication. [The factor of 
Z 2 automatically ensures that F^{k 2 = /i 2 ; P = 0) = 1.] We illustrate the realisation of the 
first two identities, Eqs. (|34]) and (|35D , in Fig. || The remaining two identities, Eqs. ( |36|) 
and Q3"7|), are realised in a similar fashion. 



V. SUMMARY AND CONCLUSION 



With renormalisation considered explicitly, we have studied the Dyson Schwinger equa- 
tion [DSE] for the dressed-quark propagator, S(p); the homogeneous, pseudoscalar me- 
son Bethe-Salpeter equation [BSE], which provides the meson's Bethe-Salpeter amplitude, 
r(fc; P), with P the total momentum of the bound state; and the inhomogeneous BSEs for the 
fully-amputated axial- vector and pseudoscalar vertices, T 5fJi (k; P) and T 5 (k; P), respectively. 
Independent of assumptions about the form of the quark-quark scattering kernel, K(q, k; P), 
we have elucidated the manner in which the axial- vector Ward-Takahashi identity [AV-WTI] 
correlates these Schwinger functions and provides important and phenomenologically useful 
constraints. 

We demonstrated that the axial- vector vertex contains a pseudoscalar, bound state pole 
contribution whose residue is the meson's leptonic decay constant, Eq. (p3|). The pseu- 
doscalar vertex also contains such a pole term but in this case the residue is related to an 
"in-meson" quark condensate, Eq. (pE2"D, which is equal to the vacuum quark condensate 
in the chiral limit, Eq. ([43]). The AV-WTI necessarily entails an identity between these 
residues, Eq. ([41]), which is valid independent of the current-quark mass of the bound state 
constituents. The expression commonly known as the Gell-Mann-Oakes-Renner relation is 
a corollary of this identity. The AV-WTI also places constraints on the form of the pseu- 
doscalar meson Bethe-Salpeter amplitude [BSA], which necessarily involves quantitatively 
important terms proportional to 757 • P and 7 • kk ■ P that have hitherto been neglected in 
phenomenological studies of spectra, and production and scattering processes in QCD. In 
the chiral limit these constraints take a very simple form: Eqs. (|3~4"D-([3~7]). 
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In Sec. [IV B 1| we illustrated these identities and constraints in numerical studies using a 



rainbow-ladder truncation for K(q, k; P) motivated by the Abelian approximation, Sees. [II B 



and [Ti rB|, and defined by an Ansatz for the dressed-quark-quark interaction, Fig.[l]. The 



model thus defined preserves the one-loop renormalisation-group behaviour of QCD, as is 
clear in our numerical solutions. This aspect of our study facilitated an elucidation of 
the dominant, ultraviolet [large- k 2 } behaviour of the scalar functions in the pseudoscalar 
meson BSA, Fig. |7]. Employing these results in a calculation of the electromagnetic pion 
form factor yields Q 2 F n (Q 2 ) = constant, up to lnQ 2 -corrections; a direct consequence of 
the model-independent result: °F H (k 2 ) ~ k 20 G K {k 2 ) ~ l/[A; 2 (ln ^/Aqcd) 1 ^™] Q. The 
pseudoscalar amplitude E n (k] P) does not contribute to the asymptotic form of F^Q 2 )^ 
In the course of these numerical studies we explored the feasibility of two methods for 
solving the homogeneous Bethe-Salpeter equation with dressed-quark propagators, deter- 
mined numerically, and without a three-dimensional reduction: (A) treating the problem 
directly as a multidimensional integral equation; or (B) employing a Chebyshev expansion 
of the scalar functions in the BSA in order to obtain a system of one-dimensional integral 
equations. We favoured (A) because it provides a simpler numerical procedure, the pro- 
gramme runs quickly and computer memory usage was not a consideration. Nevertheless, 
we observed that the Chebyshev expansion converged quickly, with at most two moments 
being necessary to reproduce the solution obtained using (A), Sec. [IV B 1| . Importantly, we 
saw that retaining the fc-P-dependence of the scalar functions in the meson BSA is all that is 
necessary to ensure that physical observables are independent of the momentum partitioning 
parameter that appears in the definition of the relative momentum, which is arbitrary in 
covariant studies. 

We have demonstrated that the Goldstone boson character of flavour- nonsinglet, pseu- 
doscalar mesons is no impediment to their description as bound states in QCD. It is only 
important that, in developing the bound state equations, the axial- vector Ward-Takahashi 
identity be preserved explicitly. This identity necessarily entails relations between the n- 
point Schwinger functions (propagators and vertices) relevant to the bound state equation; 
relations which provide for the manifestation of Goldstone's theorem. In taking account of 
this, one can in principle construct a single kernel for the Bethe-Salpeter equation that will 
provide a uniformly good, qualitative and quantitative description of the properties of all 
mesons. 



We note that the calculation of the elastic pion form factor involves two pion Bethe-Salpeter 
amplitudes, which is qualitatively different to the electroproduction studies of Ref. ||, which involve 
only one meson Bethe-Salpeter amplitude. We therefore expect that the conclusions of Ref. || 
will not be qualitatively sensitive to the omission of all vector meson Dirac amplitudes other than 
7 M since, in that calculation, the Bethe-Salpeter amplitude acts only to restrict the support of the 
integrand and the asymptotic behaviour of the cross sections is determined by the behaviour of 
the quark propagator. 
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FIG. 1. A comparison of Q(k 2 ) /k 2 in Eq. ( |23| ) obtained using the best-fit parameters of Eq. (|2 
(solid line), with 4-7ra(fc 2 )//c 2 in Eq. (||) (dashed line). The obvious infrared enhancement is qualita- 
tively and semiquantitatively in agreement with that inferred in the gluon DSE studies of Ref. |2 
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FIG. 2. The renormalised dressed-quark mass function, M(p 2 ), obtained by solving Eq. (50) 
using the parameters in Eq. (24): u/d-quark (solid line); s-quark (long-dashed line); and chiral limit 
(dashed line). The renormalisation point is fi = 19GeV. The intersection of the line M 2 (p) = p 2 
(dotted line) with each curve defines the Euclidean constituent-quark mass, M E . 
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FIG. 3. The zeroth Chebyshev moment of En(k;P): chiral limit (dotted line); 7r-meson (solid 
line); If -meson (long-dashed line); and fictitious, ss bound state (dashed line), rjp = \ in each 
case. For ease of comparison the BSAs are all rescaled so that °En(k 2 = 0) = 1. 
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FIG. 4. The zeroth Chebyshev moment of Fi{{k;P): chiral limit (dotted line); 7r-meson (solid 
line); K- meson (long-dashed line); and fictitious, ss bound state (dashed line), rjp = \ in each 
case. For ease of comparison the BSAs are all rescaled so that °En(k 2 = 0) = 1. 
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FIG. 5. The zeroth Chebyshev moment of Gn(k;P): chiral limit (dotted line); 7r-meson 
(solid line); fictitious, ss bound state (dashed line); and of Gff(fc; P) (GeV -1 ) for the if-meson 
(long-dashed line), rjp = \ in each case. For ease of comparison the BSAs are all rescaled so that 
°E H (k 2 = 0) = 1. 
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FIG. 6. The zeroth Chebyshev moment of Hn(k;P): chiral limit (dotted line); 7r-meson (solid 
line); K- meson (long-dashed line); and fictitious, ss bound state (dashed line), rjp = \ in each 
case. For ease of comparison the BSAs are all rescaled so that °En(k 2 = 0) = 1. 
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k 2 (GeV 2 ) 

FIG. 7. The asymptotic behaviour of the zeroth Chebyshev moments of the functions in 
the 7r-meson BSA: f n °ET T (k 2 ) (GeV, solid line); f-K°F n (k 2 ) (dimensionless, long-dashed line); 
k 2 ^G-nik 2 ) (dimensionless, dashed line); and k 2 f n ^^(k 2 ) (GeV, dot-dashed line). The mo- 
mentum-dependence is identical to that of the chiral-limit quark mass function, M(p 2 ), Eq. ( |59| ) 
(GeV, circles: o). For other pseudoscalar mesons the momentum dependence of these functions is 
qualitatively the same, although the normalising magnitude differs. 
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FIG. 8. An illustration of the realisation in our model of the identities Eqs. ( |34| ) and (|35|), 
which are a necessary consequence of preserving the axial vector Ward-Takahashi identity. We 
plot f°EH(k;Q) (GeV, solid line); Fn(k;Q) (dimensionless, dashed line); f°Fif(k;Q) (dimensionless, 
dotted line); and -Fft(fc;0) + 2 /#.?#(&; 0) (long-dashed line). In each curve the plotted points 
represent the right-hand-side of these equations as obtained in the solution of the chiral-limit 
quark DSE: B(k 2 ) (GeV, o); A(k 2 ) (dimensionless, o). 
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TABLES 



All 51 TY1 T\ 1 1 1 11 H £*<i 
iT.ll dllljjll U UUCD 




7T 






chiral limit 






ss 








u 


Tl H 


mo 


f° 


n H 


m sS 


fss 


n H 


Method (A) 


0.1385 


0.0924 


1.01 


0.0 


0.0898 


1.00 


0.685 


0.129 


1.00 


Uq only 


0.136 


0.0999 


0.95 


0.0 


0.0972 


0.94 


0.675 


0.137 


0.95 


?7o and U\ 


0.1385 


0.0925 


1.00 


0.0 


0.0898 


1.00 


0.685 


0.129 


1.00 


E only 








Method (A) 


0.105 


0.0667 


1.82 


0.0 


0.0649 


1.81 


0.512 


0.092 


1.68 


Uq only 


0.105 


0.0667 


1.82 


0.0 


0.0649 


1.81 


0.513 


0.092 


1.69 


£, F 








Method (A) 


0.136 


0.0992 


0.95 


0.0 


0.0965 


0.95 


0.677 


0.137 


0.95 


Uq only 


0.136 


0.0992 


0.95 


0.0 


0.0965 


0.95 


0.678 


0.138 


0.95 


£, F, G 








Method (A) 


0.140 


0.0917 


1.01 


0.0 


0.0891 


1.00 


0.688 


0.128 


1.01 


Uq only 


0.136 


0.0992 


0.95 


0.0 


0.0965 


0.95 


0.678 


0.138 


0.95 


Uq and [/"i 


0.140 


0.0917 


1.01 


0.0 


0.0891 


1.00 


0.689 


0.128 


1.01 



TABLE I. Calculated values of the properties of light, pseudoscalar mesons composed of a 
quark and antiquark of equal- mass. The mass (m° p = 0.1385) and decay constant (/| xp = 0.0924) 
are in GeV, TZh is dimensionless. With the exception of the calculations that retain only the 
zeroth Chebyshev moment, labelled by u Uq only", the results are independent of the momentum 
partitioning parameter, r\p. 
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All amplitudes 


Vp 


— u.ou 




VP 


n ot; 

— u.zo 




VP 


n nn 

— u.uu 






m K 


fK 


K K 


m K 


fK 


K K 


m K 


fK 


K K 


ivicLiiuu. yi\ j 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


1 I f\ f~\ T 1 1 ~\T 

uq oiiiy 


0.469 


0.117 


0.96 


0.482 


0.117 


0.95 


0.475 


0.113 


1.02 


TT r . nrir i rr. 
u o aiiu. cy ^ 


0.500 


0.111 


1.00 


0.497 


0.109 


1.01 


0.498 


0.110 


1.00 


U\ and t/2 


0.497 


0.109 


1.01 


0.497 


0.109 


1.01 


0.496 


0.109 


1.01 


only 








IVltJLllUU yr\. J 


0.430 


0.079 


1.55 


0.430 


0.079 


1.55 


0.429 


0.076 


1.55 


1 1 All 1 T 7" 

Co omy 


0.380 


0.077 


1.54 


0.401 


0.076 


1.51 


0.415 


0.073 


1.55 




0.439 


0.089 


1.52 


0.430 


0.078 


1.55 


0.431 


0.076 


1.57 


U , Ui and [7 2 


0.430 


0.078 


1.55 


0.430 


0.078 


1.55 


0.427 


0.076 


1.55 


F 








Method 


0.587 


0.17 


0.79 


0.557 


0.14 


0.86 


0.533 


0.11 


0.94 


T/n onlv 


0.505 


0.12 


0.82 


0.518 


0.11 


0.86 


0.512 


0.11 


0.96 


/7n and f /i 








0.556 


0.14 


0.86 


0.537 


0.12 


0.94 


Uo, U\ and U2 


0.583 


0.16 


0.79 


0.557 


0.14 


0.86 


0.532 


0.12 


0.93 


E, F, G 








Method (A) 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


U only 


0.471 


0.116 


0.96 


0.484 


0.116 


0.95 


0.477 


0.112 


1.02 


Uo and U\ 


0.504 


0.110 


1.00 


0.500 


0.108 


1.01 


0.502 


0.109 


1.00 


Uq, U\ and L^2 


0.500 


0.108 


1.01 


0.500 


0.108 


1.01 


0.499 


0.108 


1.01 



TABLE II. Calculated properties of the i^-meson for various values of the momentum par- 
titioning parameter, r\p; "— " means that no bound state solution exists in this case. The mass 
( m A" P = 0.496) and decay constant (/^ p = 0.113) are in GeV, IZk is dimensionless. 
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